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We present numerical results for gauge theory in 2d and 4d spaces involving a non- 

commutative plane. Simulations are feasible thanks to a mapping of the non-commutative plane 
onto a twisted matrix model. In d = 2 it was a long-standing issue if Wilson loops are (partially) 
invariant under area-preserving diffeomorphisms. We show that non-perturbatively this invari- 
ance breaks, including the subgroup SL(2,R). In both cases, d — 2 and d — 4, we extrapolate our 
results to the continuum and infinite volume by means of a Double Scaling Limit. In d = 4 this 
limit leads to a phase with broken translation symmetry, which is not affected by the perturbatively 
known IR instability. Therefore the photon may survive in a non-commutative world. 
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1. Non-commutative U(l) gauge theory 

The positions in a non-commutative (NC) Euclidean plane correspond to the spectra of Her- 
mitian operators jq, #2, with a non- vanishing commutator 

[xn,x v ] = j0 MV = iQfyv ■ (1-1) 



We treat the non-commutativity parameter 6 as a constant. Relation ( |1.1| ) implies a spatial uncer- 
tainty of the form Axi Ax2 ~ 0, which can be interpreted as the event horizon of a strong gravitation 
centre. In fact, if it has been argued that attempts to merge quantum theory with gravity lead quite 
generally to such a spatial uncertainty [|I|], which corresponds to a NC geometry. 

In quantum field theory the spatial uncertainty gives rise to non-locality over a range of 
0(Vd). A related consequence is the notorious "UV/IR mixing" []2|]: nested singularities can 
be UV divergent in one momentum component and IR divergent in another component p v . Due 
to this property the perturbative treatment is extremely involved. 

Hence it is strongly motivated to take a fully non-perturbative approach. As in commutative 
field theory it relies on the lattice regularisation (for a review, see Ref. [[§]). Although we do not 
have sharp points as lattice sites, a (fuzzy) lattice structure can be imposed by the operator identity 

exp I i — ) = 1 , (1.2) 



a 

where a is the lattice spacing. Along with the usual periodicity of the momenta over the Brillouin 



zone, identity (L2) entails ^Op^i 6 ^ • For fixed parameters a and 6 we infer that the lattice is 



automatically periodic, in striking contrast to the commutative lattice. 
On a periodic N xN lattice one readily identifies 

B = -No 1 . (1.3) 

Therefore we extrapolate to the Double Scaling Limit (DSL) 

{ a^O and iV^» } at No 2 = const. (1.4) 

The DSL leads to a continuous NC plane of infinite extent. The requirement to take the UV and IR 
limits simultaneously in a controlled manner is again related to the UV/IR mixing. 
We can return to ordinary coordinates x^ if we multiply all fields by star products, 

0(x)*yA(x) := 0(x)exp ^^©^y^v) W(x) , (1.5) 

which encode the non-locality. The star commutator [xu,x v ]* := x^ *x v —x v -kx VL = i&n V suggests 
that this transition is sensible — it can be justified e.g. with a plane wave decomposition. 
In this framework we formulate U(l) gauge theory on a NC plane as 



S[A] = ^J d 2 x F^y-kF^y , Fuv = d^A v -d v A^ +ig[A fl ,A v }± . 



(1.6) 



Note that even the U (1) gauge field picks up a self-interaction term. This action is invariant under 
star gauge transformations. However, even on the lattice its direct simulation is hardly feasible (for 
instance the compact formulation would require star unitary link variables). 
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We arrive at a numerically tractable formulation based on the equivalence of this model with 
the twisted Eguchi-Kawai (TEK) model [§]. This matrix model is defined on one point with the 
action 

SrEK[U] = -NPY t Z liV 'rr[U ll U v Ulul] , j3 = l/g 2 . (1.7) 

The are unitary N xN matrices, which contain the degrees of freedom of the lattice gauge field. 
The twist factor Z21 = Z* 2 = exp(2nin/N) makes the difference from the original Eguchi-Kawai 
model — in d = 4 it avoids the spontaneous breaking of the centre symmetry at weak coupling, n 
is an integer, which we set n = (N +l)/2, so we always deal with odd matrix/lattice sizes N. 

Then the TEK can be identified with the NC U(l) lattice gauge model since the algebras are 
identical — this has been shown in the large N limit |Q] and also at finite N [|6p. Clearly, the TEK 
is numerically tractable, so that the simulations can start if we also formulate suitable observables. 
In the matrix model framework it is obvious to write down the analogue of a Wilson loop, 

WnvQ x J) ■= ^Z^Trpftiul'ul 7 } . (1.8) 

Mapping this quantity back to the lattice yields indeed the star gauge invariant term, which is 
considered the NC Wilson loop |0|. This Wilson loop is complex in general, W^ v £ (C, but the 
action is real since both orientations are summed over (and W^ v = Wy„). Hence simulations are 
possible without running into a sign problem. Such a study, including extrapolations to the DSL, 
was first presented in Ref. The planar limit coincides with U(N — > °°) lattice field theory, 
where the (real) Wilson loop follows an exact area law [g]. Although this is not the limit that we 
are interested in, it can be used to set the scale as a 2 = — ln[l — 1/ (4/3)] (for /3 > 1 /2), so that the 
DSL has to be taken roughly at a fixed ratio N/ft. 



2. Wilson loops in d = 2 : area-preserving diffeomorphisms (APDs) 

On the commutative plane, pure U (n) gauge theories are analytically soluble with geometric 
methods. Due to APD invariance, the expectation values of Wilson loops only depend on the 



oriented area [1C]. 

Contrary to original expectations, it turned out that this symmetry does not hold on the NC 
plane. In particular, perturbation theory to 0(g 4 ) and 0(d~ 2 ) revealed that it breaks down to 
SL(2,R) [11], in agreement with other considerations [|T2|]. By means of numerical simulations we 
investigated the non-perturbative extent of this APD symmetry breaking, as well as the viability 



of the subgroup SL(2,R) QX 3[] . To this end, we considered four types of non-intersecting Wilson 
loops with polygonal boundaries, generalising the form (|L8|). Prototypes are depicted in Figure [[]. 
Figure || shows results for a set of Wilson loops at N = 125, j8 = 3.91 (on the left) and N = 155, 
ft = 4.82 (on the right). In both cases the non-commutativity parameter amounts to 6 = 2.63, hence 
we see two snapshots on the way to the DSL: the right-hand-side corresponds to a larger volume 
with a finer lattice. For a fixed (dimensional) loop area the results are very similar, hence we have 
apparently reached the asymptotic DSL behaviour. At small areas we observe agreement with the 
Gross-Witten area law [^] (and therefore with the planar limit) for all shapes. On the other hand, 
at larger areas \(W) \ does not decay further and the values for the different shapes drift apart. This 
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5 



Figure 1: Examples for squares, L-shapes, rectangles and stairs, the four types of Wilson loops that we 
measured at the same area to investigate the behaviour under APD transformation, cf. Figure 2. 

shows the extent of APD symmetry breaking, which seems to persist in the DSL to the continuum 
and infinite volume. The rectangles (including squares) are related by the APD subgroup SL(2,R). 
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Figure 2: Absolute Wilson loop values of different shapes for two volumes and lattice spacing at 9 = 2.63. 
The good agreement between the two plots indicates that the APD symmetry breaking persists in the DSL. 

In fact, SL(2,R) has a higher viability as an approximate symmetry subgroup — at least up to 
moderate deformations — but it also breaks on the non-perturbative level. 



3. The fate of the photon in a non-commutative world 

We consider again a NC plane with [xi,X2] = id = const., but now we add a commutative plane 
(x3,Xt), which includes the Euclidean time. We are particularly interested in a possible 6-distorted 
dispersion relation of the photon, which could in principle be experimentally measurable. A one 



loop calculation suggests the form [ 14 ] 

PuPn 



where C = const. 



P\i ■= ®iivPv 



(3.1) 



Corresponding phenomenological data have been analysed, in particular in view of the time of 
flight of cosmic photons in Gamma Ray Burst (GRBs). In a GRB photons over a range of about 
10 s . . . 10 eV are emitted within a few seconds or minutes from a small source. In particular, eval- 
uating the times of photon arrival for 35 GRBs against the dispersion ansatz E = \p \ + E/M (where 
M is a large mass due to some "quantum gravity foam"), Ref. [ jl5| ] concluded M > 0.001Mpi anc k. 



It has been proposed to establish a bound on [|©[| by similar considerations [16], which should ad 



dress the IR singularity in eq. (3.1). However, the constant C turned out to be negative to one loop 
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Figure 3: Results for Wilson loop in NC QED4 at N /pi 2 = 20 in the commutative plane, the mixed plane 
and the NC plane (from left to right). In the first two planes (W) is real, whereas the oscillation of Re (W) 
comes along with a rotating complex phase (similar to NC QED2 /@7J. 
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Figure 4: The expectation values of open Polyakov lines of length 2a (on the left) and 4a (on the right) for 
N = 15, 25 and 35 (curves from left to right, with hysteresis at the transition to the weak coupling phase). 
At intermediate gauge coupling we recognise a phase of broken translation symmetry. 



[JlVp , which suggested that NC QED is IR unstable and thus ill-defined. In Ref. [ ]18[ ] we revisited 
NC QED non-perturbatively: we discretised the commutative plane with an L x L lattice and the 
NC plane with a TEK model of matrix size N m L. A physical scale was identified by matching 
Wilson loop expectation values (and further observable) at different j8 values. To a good approx- 
imation this led tofl« 1//3, so that the DSL is taken at fixed d <=c N/f5 2 (unlike NC QED 2 ), see 
Figure BJ. 1 

As an order parameter for translation symmetry in the NC plane we measured the open Polyakov 
line, which is *-gauge invariant and which carries momentum p, 



P fl (n) = 0>expjig 



& : path ordering, =najl : length . (3.2) 



Figure |] shows the (expected) symmetric phases at strong and at weak coupling, but a broken phase 
above j8 ~ 0.35. The upper end of the broken phase rises roughly <=c N 2 , and its hysteresis behaviour 
implies that the corresponding phase transition is of first order. 2 

This leads to the phase diagram ||, where we also mark the trajectories for DSLs with different 
values of 8. These DSL curves always lead to the broken phase, where we observe IR stability 



of all observables measured Q18|]. The perturbative result describes correctly the weak coupling 
phase, as the dispersion relation (in the commutative plane) in Figure |6| on the left shows. The 

'The precise fine tuning results, with slight deviation from this rule, are given in Ref. pl8[], 

2 The broken phase of the TEK model was mentioned earlier Jl9[ ] and later confirmed in Ref. [po|]. 
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Figure 5: The phase diagram for QED4 with a NC plane: between the strong coupling phase (fi < 0.35 j and 
the weak coupling phase we find a phase of broken translation invariance. Since the weak/broken transition 
rises like /3 °= A^ 2 , whereas the DSL for a fixed NC parameter 6 follows /3 « \^N, the DSL always ends up in 
the broken phase. 
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Figure 6: Dispersion relations E(p) for the photon in a NC space in the weak coupling phase (on the left, 
p = p\ ), and in the broken phase ( on the right, p — pj,). The former result is consistent with the perturbatively 
predicted IR instability ("tachyonic" behaviour, eq. (3.1) with C ~ 0.1285j. In the physical phase the photon 
is massless again, corresponding to a Nambu-Goldstone boson of the broken translation symmetry. 

plot on the right refers to the broken phase, which captures the physically relevant DSL. Here the 
photon dispersion is linear, and the photon can be identified with the Nambu-Goldstone boson of 
spontaneous translation symmetry breaking. 3 



4. Conclusions 

We simulated QED2 and QED4 on spaces containing a NC plane. In both cases we observed a 
stable behaviour in the DSL to a continuous NC space of infinite extent, which suggests renormal- 
isability. UV/IR mixing is manifest as a non-perturbative effect, in agreement with numerical DSL 
results for the 3d NC A0 4 model [|23|]. 

3 This interpretation is known for instance from Ref. The ansatz taken there incorporates a NC space where the 
star product is truncated in O(||0||) [122]. In that case, however, locality is restored, so the theory is altered qualitatively. 
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• In d = 2 the APD invariance of Wilson loops breaks, without any residual subgroup. Hence 
there is hardly hope for an analytic solution, but we may hope for a rich structure to be explored 
numerically [||, 

• In d = 4 the DSL leads to a phase of intermediate coupling strength and broken translation 
symmetry. This physical phase appears to be IR stable — in contrast to the weak coupling phase. 
Therefore the NC space may accommodate photons after all [18]. 
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